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Abstract 

We propose to use a Nomarski imaging interferometer to measure 
the out-of-plane displacement field of MEMS. It is shown that the 
measured optical phase arises both from height and slope gradients. 
Using four integrating buckets a more efficient approach to unwrap 
the measured phase is presented, thus making the method well suited 
for highly curved objects. Slope and height effects are then decou- 
pled by expanding the displacement field on a functions basis, and 
the inverse transformation is applied to get a displacement field from 
a measure of the optical phase map change with a mechanical load- 
ing. A measurement reproducibility of about 10 pm is achieved, and 
typical results are shown on a microcantilever under thermal actua- 
tion, thereby proving the ability of such a set-up to provide a reliable 
full- field kinematic measurement without surface modification. 

1 Introduction 

The increasing interest for micro-electro-mechanical systems (MEMS), es- 
pecially when they are used as micro- mechanical sensors |T], leads one to 
focus on the mechanical behavior of micro-objects. First, the standardized 
mechanical tests at the macro scale have been adapted to the micro one, as- 
suming an homogeneous stress or strain state. Bending [2] and tensile tests 
[3 EI, as well as fatigue or creep tests [S] are performed since years, providing 
a global kinematic response of the tested object. 

However, classical photolithography processes use visible light to transfer 
a mask onto a wafer surface. Then, light diffraction limits the achievable 
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accuracy of the geometric shape of the resulting micro-objects, and the di- 
mensions margins tends to increase compared to the dimensions themselves 
as the object's size decreases. The homogeneous stress (or strain) assump- 
tion usually satisfied when performing mechanical tests at the macro-scale 
is no longer reasonable, and one has to deal with heterogeneous mechanical 
tests. As a consequence, one has to perform a spatially resolved kinematic 
measurement instead of a global one. 

Moreover, as their size decreases, the surface on volume ratio significantly 
increases, and the behavior of micro-objects tends to be dominated by surface 
effects instead of volume ones. As a consequence, measuring the displacement 
field should avoid any surface modification or any contact. In addition, the 
surface roughness of MEMS is usually very low, so that measurement tech- 
niques involving surface-generated speckle jH] may be difficult to implement 

® 

An optical interferometric imaging set-up is then well suited to measure 
displacement fields. The polarization interferometer [S] proposed herein is 
derived from the one initially proposed by Nomarski Dealing with surface 
topographies, it has been previously used to determine the mean slope of 
tilted samples [TU] or to get an image of the roughness of polished surfaces 
using a multichannel Nomarski microscope jTH . 

After recalling the relation between the optical phase induced by the 
sample- Wollaston prism group and the measured intensity, the twofold origin 
of the optical phase is detailed. A shot-noise limited detection is described, 
and the inversion method developed to convert an optical phase change into 
an out-of-plane displacement is presented. An example is finally provided 
on the measurement of the displacement field of a thermally loaded micro- 
cantilever. 

2 Measuring a differential topography 
2.1 Experimental set-up 

The basic interferential microscopy imaging set-up is shown in Fig. ^ A 
spatially incoherent light source (LED, A = 760 nm) illuminates a polariza- 
tion beam-splitter. The beam reflected by the beam-splitter goes through a 
polarization modulator and is initially polarized at 45 ° of the axes of a Wol- 
laston prism. This splits the beam into two orthogonally polarized beams at 
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a small angle between each other. These beams are focused upon the sample 
by an objective lens. After reflection and recombination by the Wollaston 
prism, the beam goes through the polarization modulator and the polariza- 
tion beam-splitter. The transmitted beam is finally focused on a CCD array. 
The polarization beam-splitter behaves as crossed linear polarizers mounted 
at 45 ° of the axes of the Wollaston prism and of the polarization modulator. 
One finally gets the interference between two images of the sample, shifted 
by the Wollaston prism of the distance d. The resulting interference pattern 
is recorded on a CCD array (DALSA CA-D1, 256 x 256 pixels, 8 bits). 

2.2 Interference pattern obtained with a Nomarski imag- 
ing interferometer 

Let us denote p (FigJTJ the polarisation direction of the beam incident on 
the Wollaston prism. The Wollaston prism shear-direction is denoted y. The 
vector a denotes the polarisation direction of the beam impinging on the CCD 
array. The orthogonal directions, in the prism's plane, are denoted q, x, b, 
respectively. One uses a light emitting diode, so that we denote E the am- 
plitude of the non polarized wave impinging on the polarizing beam-splitter. 
If t p is its amplitude transmission factor, and e p the attenuation factor for 
the (ideally) suppressed component, the Wollaston prism is illuminated by 
two orthogonally polarized beams, which electric fields are E pp p and E pq q, 
with 

E pp = E t p (1) 

E pg = E € p tp (2) 

For a non polarized light source (LED), these two beams are incoherent, 
and should be treated separately. e p t p is the transmission factor in the stop 
direction, so that e p = is the perfect polarizer case. For each beam, the 
Wollaston prism splits the beam into two orthogonally polarized beams (i.e., 
a x and a y component), and the light goes through the path (objective 
- sample - objective - Wollaston prism). Let us consider that this results 
in a phase shift <fi between the x and the y components. Let us denote t n 
the transmission factor of the path (Wollaston prism - objective - sample - 
objective - Wollaston prism), t a the analyser's transmission factor, and e a t a 
its further attenuation factor in the stop direction. The total intensity / 
impinging locally on the CCD array is the sum of the intensity arising from 
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the two incoherent beams 




(3) 



Assuming that (x, p) = (: 
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The contrast of the interference pattern 
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equals 1 for perfect polarizers (e a = e p = 0), and decreases when e a or 
e p increase. Adding an error on the orientation of the polarization beam- 
splitter leads to the same expression, so that the equation (jlj) is considered 
general enough to describe real interference patterns. A typical interference 
pattern obtained in water with two 70 micrometers long micro-cantilevers 
and a shear distance d ~ 50 is shown on figure |21 The two sheared 
images are clearly distinguishable. The optical phase range covers almost 15 
interference fringes. The closely packed interference fringes, as well as the 
quite short correlation length of the used light source (almost 15/wn) reduce 
the contrast of the interference pattern, thereby limiting the quantitative use 
of obtained phase map. Let us assume that the optical phase introduced 
by the path (Wollaston prism - objective - sample - objective - Wollaston 
prism) may be decomposed in a term O arising from the Wollaston prism 
and a contribution m arising from the object 



The section 12751 exhibits the phase directly arising from the Wollaston prism, 
whereas the section E3] exhibits the phase arising from the topography of the 
sample. 

2.3 Optical phase arising from the Wollaston prism 
2.3.1 Optical path functions 

Let us consider for simplicity an "ideal" Wollaston prism, which geometry is 
described in the plane defined by both the optical axis of the system and the 



= O + 0, 
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y direction on figure 01 Fig. |3] shows the decomposition of a ray impinging 
orthogonally on the prism into two emerging rays and the angles definition. 
Snell-Descartes' laws at the interface between the two half-prisms read 



n E sm(f3™) = n sm(6) (9) 
n sm(p TE ) = n E sm(9) (10) 

where n is the ordinary refractive index of the used birefringent material, and 
n E is the extraordinary one. Assuming an ambiant media with a refractive 
index equals to 1, Snell-Descartes' laws for the exit face of the prism read 

sin(/3™) = n E sm((3™-6) (11) 
sin(^f) = n sm{(3 TE -6) (12) 



If all the angles are small 



02?-0%^2{n o -n s )9 (13) 



The emerging rays appear to split on a plane, called plane of apparent split- 
ting (PAS), which is here located inside the prism. The use of a modified 
Wollaston prism ^3] allows one to move the PAS outside the prism. This 
plane is considered to be perpendicular to the figure's plane. Its position is 
one of the fundamental characteristics of the prism. Investigating the imag- 
ing properties of the system, let us denote a.\ the angle of an impinging ray 
with respect to the normal of the entrance interface in the figure's plane. 
Considering any dependance on a% is then moving in the field of view along 
the shear direction. The optical path travelled by the TE (resp. TM) polar- 
ized ray through the prism l TE (y,ai) (resp. l™(y,ai)) when the apparent 
splitting occurs at position y on the PAS (the origin will be defined later) 
depends on a±. Assuming that all the angles are small, 

These two functions will be referred as optical path functions. 
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2.3.2 Optical phase when the PAS matches the objective rear 
focal plane 

Assuming that the PAS of the prism matches the objective rear focal plane, 
Fig. EJpresents the ray tracing for the two emerging rays of Fig. El These two 
emerging rays intersecting the PAS at the y pas cross the PAS, after reflection 
on a plane object, at the the same point Ur^te = Vr,tm when the object is 
orthogonal to the lens axis. This last point is symmetric of the first with 
respect to the lens axis. Let us thus define the median plane of the prism ir m 
as orthogonal to both the figure's plane and the PAS, passing by the O point 
(see Fig. |SJ). This point is defined to satisfy the condition EO = OS. Let 
us then consider that the distance between the 7i m plane and the symmetry 
axis of the lens is T w . Any point of the PAS may be described either by 
the abscissa y in the prism's frame (the origin is at the intersection with 
ix m ) or by the abscissa y = y — T w (relative to the symmetry axis of the 
lens). In the prism's frame, the rays emerge from the PAS at y pas = T w — y t , 
and the reflected rays go through the PAS at Vr^te = Ur,tm = T w + y t . The 
optical paths Zj r E and /™ travelled to the sample and back are deduced from 
l TE (y, ai) and l™(y, ati) and allow one to compute the optical phase 

<hro = ^(C-C M ) = ^ U*b - n )T w + harm + ^(^ - 1) - lY) 

(15) 

The first term in ()15j) no longer depends on y t but on the position of the prism 
with respect to the symmetry axis of the lens T w , and is homogeneous in the 
field of view. If the birefringence of the used material is denoted % — n = 
n Q e, the additional optical phase difference (second term) proportional to the 
incidence angle a% is shown to grow as e 2 whereas the first term scales as 
e, so that the added optical phase (fiwo ma Y be considered homogeneous in 
the field of view. For a quartz-made Wollaston, e ~ so that the added 
optical phase may be considered homogeneous in the field of view. This 
is no longer true when the PAS doesn't match the rear focal plane of the 
objective, because the reflected rays no longer cross the PAS at the same 
point (see Fig. HJ. Assuming that the PAS remains parallel to the rear focal 
plane, at a distance <5pAs and using the same method, the additional optical 
phase difference 4>pas is 

87T 

(Ppas - ^OatSpAsiriE - n Q ) (16) 
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thus introducing a linear phase along the Wollaston shear direction. 



2.4 Optical phase arising from the object topography 

2.4.1 Optical phase arising from height variations 

The previous ray tracing is shown in the case of a stepped sample (height 
Az) in Fig. Thanks to the Fermat's principle, tilting the sample doesn't 
induce an extra phase shift in the objective-sample path (i.e., regardless of 
the phase shift induced by the Wollaston prism). The dot line is for the 
ray reflected under the previous conditions. The optical phase difference 4>h 
arising from the step corresponds to a travel along [AB] and [BC], where 
C is the orthogonal projection of A on [-B-D]. In a medium which refractive 
index is n, this optical phase difference reads 

271 4:71 

<p h = —n x (AB + BC) = —nAz cos(a) (17) 
A A 

where a is the incidence angle on the object. If the numerical aperture of the 
objective is low enough, a remains small and the above expression expands 

AnnAz 

4>h - — t — (18) 

As the numerical aperture of the objective increases, the previous expansion 
is no longer valid and the equation (jTSjl is replaced, for numerical apertures 
less than 0.3, by 

4irnAz 

* = — (19) 

where i is a scale parameter depending on the numerical aperture (see ap- 
pendix EJ. 



2.4.2 Optical phase arising from the local surface orientation 

In section D2.4.11 the two orthogonally polarized rays cross on the PAS. This 
is no longer true if the two rays experience a different surface orientation. 
Let us denote by r y TE (resp. 7™) the local surface orientation experienced 
by the TE (resp. TM) ray. For a ray emerging from the PAS at the position 
y (in the prism's frame), the reflected rays cross the PAS at the abscissa 
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(figure IHJ) 



_ -y + /tan(2 7 rg ) 

" l + ftan(2 7 ^) (20) 

-y + /tan(2 7 ™) 
y *™ " l + ftan(2 7 ™) ( } 

Where / is the focal length of the objective lens. The optical path travelled 
by the two rays reads (with a\ = 0) 

l T ar E ~ fZ(T w + y,0)+l TE (T w + y R , TE ,0) (22) 
I™ ~ /™(T w + 2/ ,0) + F M (T w + ^ TM ,0) (23) 

and the optical phase difference reads (to the first order with respect to the 
surface orientation) 

<Pw = y = y(n E -n o )tan(0) ^1 + (j)' I /(7t E +7tm) (24) 

One should highlight that this optical phase difference doesn't arise from the 
sample itself, but from the optical phase imbalance it introduces when the 
rays travel back through the Wollaston prism. Moreover, this phase difference 
depends on the abscissa y, and then on the angle a. It's also worth noting 
that a global tilt of the specimen (i.e., ^te = Itm ^ 0) induces an extra 
optical phase. It arises from the previous remarks that the relation between 
the local surface orientations and the induced optical phase depends on the 
prism-objective group. This relation may be rewritten 

4>w = 4>ori + (fruit = — [iTE ~ Itm) + 4>uit (25) 

(77 



with 



Hilt 



<9 7 



Itm (26) 



The scalar then describes the prism-objective group. It may be obtained 
when measuring the optical phase when tilting a reasonably flat mirror. The 
figure [7| shows the result of such a calibration for a quartz Wollaston prism 
[9 ~ 18°) and a water immersion objective lens (focal length 18mm, NA 0.3, 
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used to obtain the interference pattern in FigJ2J). The relation is linear, and 
the fitted slope is 

- 9<Pw - 1.1 x 10 3 (27) 



^7 

This simple experiment allows one to obtain the sensibility to a local 
orientation gap by tilting the whole sample. 

Finally, the total measured optical phase difference <fi may be written as 
the following sum 

= 4>wo(T w ) + (j)pAs(dpAs) + 4>h(Az) + (ftwilTE, 7tm) (28) 

where the two first terms (related to 0o) account for the path difference 
introduced by the Wollaston prism itself : an homogeneous term <fiy/o and 
a linear term along the shear direction <f)pAs which vanishes when the PAS 
matches the rear focal plane of the objective lens. The two last field terms 
(related to <f> m ) originate from the surface topography : height variations 0/, 
and slope field, <fiw including the contribution of the mean tilt of the sample 
4>tm as well as slope variations or j. 



3 Phase integration, retrieval and unwrap- 
ping 

3.1 Principle 

The optical flux collected by the pixel indexed by (l,m) on the CCD matrix 
can be formally written according to (@J) as 

I(l,m,t)=I + AcoB[<f>(l,m)+'if) mod (t)], (29) 

where <f)(l, m), which is to be determined, is the optical phase introduced by 
the sample and the Wollaston prism. The phase modulation introduced by 
the photoelastic modulator reads 

tpmod{t) = tpo sin(27r f mod t + 9 mod ). (30) 

The angles ipo et Q m od are two parameters that can be chosen among many 
couples. The algorithm to obtain 6 uses four integrating buckets |TH] : ifT = l/f mo d 
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is the modulation period, four images of the interference pattern can be cap- 
tured during the period T, so that each image results from the integration 
of the optical flux during a quarter of one period. One obtains four images 
E p , for p = 1,2,3,4 

E v = I(t)dt (31) 

J (P-1)T 
4 

Combining with Eq. ()29jl allows one to obtain 
T 

E p = -(I Q + A/ (V>) cos(0)) 

TAcos(0) ^ J 2 „(^o) r . / LO fl x . , , i\ , n n \ 
+ > — [sm(npn + 2nO mod ) - sm{nvp - 1)tt + 2nO mod ) 

n=l 

TAsin((f)) J 2n +i(^o) 



7T ^ 2n + 1 

?1=0 



(cos(^(2n + l)(p - 1) + (2ri + 1) 



'mod ) 



-cos(|(2n + l)p+(2n + l)^ w )) (32) 

where J n is the first kind Bessel function of n order. The images E p depend 
on J , cos(0) and sin(0), so that using four independent images provides 
enough information to recover <p. Classical algorithms ^3] use particular 
linear combinations 



S s = -(E 1 -E 2 -E 3 + E i ) = ^T sS m(<j ) ) 
S c = -(E 1 -E 2 + E 3 -E,) = ^r c cos((f>) 



(33) 



with 



r s = E~o(-l) n %# 1 sin[(2n + l)^ 



modi 



r c = E^=o^# 1 ^[2(2n + l)& 



modi 



(34) 



The optical phase is then usually recovered by using an "atan2" function 
with the arguments provided by the equations (J33j) with a (if) , 6 mod ) couple 
satisfying Y s = T c . The indicator 

X 2 = £ 2 + £ 2 (35) 

is then independent of the phase <fi and defines the intensity of the signal 
experiencing the phase <fi. As a consequence, the phase obtained from the 



10 



previous algorithm is reliable provided that T 2 is high enough, 



i.e.. 



if the 



fringe contrast is sufficient. This condition may be difficult to satisfy when 
the topography is described by a large number of closely packed interference 
fringes (see Fig|2J). It's worth noting that the set of four images provides an 
over-determinated set of data. This redundancy is then exploited to provide 
a more reliable phase measurement. 



3.2 Least-square phase retrieval 

For each pixel of the CCD array, the set of Eq. 
linear system of equations 

MP = E 



may we rewritten as a 
(36) 



where the parameters vector P reads ((•)* denotes the transpose of (•)) 



TI TA TA . 
— , cos(0), sm(0) 

4 7T 7T 



and the images vector 



E = [Ei, Ei, E 3 , £4] 
The matrix M is built from the modulation parameters 



(37) 
(38) 



M 



1 c(l,ll) ,9mod) s(l,lf) ,9mad) 

1 c(2, -00, Qmod) s(2,ij ,9 mod ) 

1 c(3, 00, Qmod) s(3,tp ,9 mod ) 

1 c(4, O , Omod) s(A,^ ,6 mod ) 



(39) 



with 

c(p,ij ,9 mod ) 

S(P, 00, Qmod) 



n=l 



2n 



sm(n(p — 1)tt + 



E 

n=0 



J2n+l(^o) 

2n + 1 



,7T 



x ( cos(-(2n + l)(p - 1) + (2n + 1)^) 



,7T 



cos(-(2n + l)p+(2n + l)^ mod ; 



(41) 
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The matrix M is then independent of the considered point. For each pixel 
of the CCD array, the four images describe the vector E, and the solution 
parameters vector P so i is obtained as a minimizer of 

^(P) = (MP - E)* {MP - E) (42) 

and is the solution of the square linear system 

M'MPsoi = M'E (43) 

The couple ^ TAco ^ s °^ , rAsi °(^°') j j s then extracted, and used as the argu- 
ment of a standard "atan2" function to provide a less corrupted value of the 
phase. The figure |H] shows the phase map obtained from the scene presented 
in Fig. |21 using the least-square algorithm. 

3.3 Phase unwrapping 

The use of the phase integration technique also allows one to reconsider the 
phase unwrapping problem. In the previous section, the information derived 
from the experiments is a couple (X,Y) proportional to (cos(0), sin(0)) 

(X,y) = C(cos(0),sin(0)) (44) 

where C = ^-^T S)C when the linear combinations (J3~3*)) are used, and C = ™ 
when the least-square algorithm is used. In the first case, the couples (X, Y) 
and (E c , S s ) are equal. For a sake of simplicity, let us consider that a wrapped 
value of the phase (ft p is then obtained by using 

(ftp = atan2(Y, X) (45) 

(ftp lies then in [— 7r,7r]. Recovering the unwrapped value of the phase may 
turn into a brain-racking task when the phase-map is poorly discretized and 
the phase jumps are densely packed. Numerous algorithms are available in 
the literature, based either on local phase unwrapping ^B] or on global phase 
unwrapping [T7j. Contrary to most of these algorithms, which use a wrapped 
phase map as an input, we propose here to use the two fields X and Y. 

The figure |H1 shows a typical case of phase jump between two adjacent 
points 1 and 2. Their true phase <p\ and 02 is under scrutiny. As X\ ~ X2 7^ 
and Y{Y 2 < 0, using the "atan2" function induces a phase jump when moving 
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from 1 (<p p i — 7r) to 2 (0 P)2 — — 7r). Let us then define the frame 

obtained by rotating the (X, K) frame until the point 1 (X±, Y\) belongs to 

the (Q, Xb) axis. The phase 02 — 0i of the point 2 in the frame (Xb, Yb) is 
then obtained according to 

, r , , ! -X 2 sin(<l> 1 ) + Y 2 cos((p 1 ) 

tan 2 - 0i = — (46) 

A 2 cos(0i) + Y 2 sin{(pi) 

providing the gap between the true phase values 0i and 02, thus defining 
(if 1 and 2 are for adjacent pixels), the true phase gradient, modulo 2n. 
The couple (X 2 ,Y 2 ) is obtained of the four images using one of the two 
algorithms previously described. cos(0i) and sin(4>x) are deduced from the 
wrapped value 4>p,i- The phase gradient may then be integrated to provide 
a true phase map. The figure E3 shows the phase map obtained from the 
wrapped phase map shown in FigJHl using the described phase unwrapping 
technique. 



3.4 Reproducibility of the phase measurement 

The phase measurement reproducibility is assessed by measuring twice the 
phase map arising from the same differential topography of a reflective object. 
One gets two phase fields 0_(/, m) and 0+(/, m). Assuming that each of these 
fields is the sum of deterministic part 0d(/,m) and of a random part 

<f>-(l,m) = (j) d (l,m) + b_(l,m) (47) 
0+(/,m) = <j> d (l,m) + b+(l,m) (48) 

the difference between the two phase fields provides then a realization of the 
difference between two realizations of the noise b(l,m) 

(0 + -0_)(/,m) = (6 + -6_)(/,m) (49) 

The figure ^2 shows a typical probability density of the variable (0 + — 
0_) obtained with the set-up described in Section 12.11 The dots are the 
experimental distribution, whereas the solid line is a least-square fit using a 
Gaussian, zero- mean, distribution. The agreement is excellent, and allow us 
to consider (0 + — 0_) as a random real number of variance 2a 2 , if b is a random 
variable which variance is a 2 . The figure IT2*1 shows the evolution of \/2a 2 
(converted to heights variation assuming (fTHjl ) as a function of the exposure 
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time t exp used to form each intensity image. The experimental variance 

_ i 

grows as t ex 2 p , until the exposure time reaches several tens of seconds. The 
reproducibility is therefore controlled by the exposure time, and the achieved 
level is (converted to height variations) close to 10 pm. 

4 Recovering the displacement field 
4.1 Principle 

The calculations presented in sections 22.31 and 22.41 allow one to derive an 
expression for the measured phase as a function of the surface topography. 
Assuming that the pixel size (in the object plane) is small enough comparing 
to the object's size, one proposes for a sake of simplicity to formulate the 
inversion problem in the object's plane. Inserting Eq. (JT9jl and Eq. (|2l)j) into 
Eq. (J25j) allows one to write the total measured phase difference at the point 

<p(x,y) = <p (x,y)+<f> 1ia t + ^(,x,y) (50) 

Assuming that the shear direction is parallel to the y axis, and that the 
local surface orientation 7 in Eq. (}25|) is given by the first derivative of 
the topography, the relation between the measured information $ and the 
topography reads, for a shear distance d 

~, . 4™ / d d.\ d(pw ( dz d s dz . d \ 

V) = [z{x, y + -)- z(x, y--)j [^(x, y+ -)--( x , y --)) 

(51) 

where z(x,y) is the topography. To discriminate between the phase aris- 
ing from the slope variations and those arising from height variations, one 
proposes to expand the topography on a functions basis 

z (x, y) = ^2 VsZs(x, y) (52) 

s 

so that the reference (initial) measured phase field reads 

(f) r ef{x, y) = (f>o{x, y) + (ptut + ^2 P&'fa y) ( 53 ) 

s 

One should highlight that the shearing amount d is here introduced explicitly 
to compute the functions basis <& s (x,y) by inserting Eq. (J52j) into Eq. (J51|) . If 
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the topography is subjected to an out-of-plane displacement field w(x), 4> re f 
is changed into (f> w and the new topography is described by 

z(x, y) + w(x, y) = ^(/x, + v,)z,(x, y) (54) 

s 

so that the displacement field is also expanded on the same functions basis 

w ( x , y) = ^2 u s z s( x i v) ( 55 ) 

s 

The knowledge of the numerical coefficients involved in the definition (|51|) of 
<&(x,y) combined with the choice of a functions basis allows one to recover 
the displacement field by computing the v s coefficients minimizing 

vM = J (y2v s $ s (x,y) - {(j) w - (f) ref )(x,y) \ dxdy (56) 
that is solving the linear system 

Mv = F (57) 



with 



N si = j <$> s (x,y)<$>i(x,y)dxdy (58) 
$ s (x,y)((f) w - (j) ref )dxdy (59) 



4.2 Example 

The micro-cantilevers shown on FigEJ are placed into a fluid cell, filled with 
milliQ water which temperature is controlled thanks to a feedback loop con- 
trolled Peltier device. These cantilevers are multi-layer cantilevers, made 
of a silica layer (770 nm), a titanium layer (20 nm) and a gold layer (50 
nm). These cantilevers are then subjected to a bimaterial effect, since the 
coefficient of thermal expansion of the gold layer is almost ten the one of 
the silica layer. The optical set-up is also used with the Wollaston prism- 
objective couple calibrated in Section 12.4.21 A reference phase map of the 
cantilever is captured at 22.6 ° C. The temperature is then increased to reach 
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24.1 ° C, and a new phase map is captured. The phase gradients are com- 
puted according to the algorithm described in Section 33.31 Recovering the 
two-dimensional phase map from the phase gradients is an over-constrained 
problem, so that it is possible to avoid some unreliable data. Reliable data 
are then located where the indicator defined by Eq. (j33j) is greater than a 
user-defined threshold. The figure IT31 shows the measured phase map change 
when the cantilever is submitted to bimaterial effect. One distinguish the 
substrate, which is not subjected to any modification. One then remark the 
phase change is not homogeneous across the cantilever. 

The displacement field is then recovered using the algorithm described in 
Section 44.11 The function basis is chosen to be able to represent the expected 
mechanical effects. As a consequence, one chooses cubic hermite polynomials 
[T8] by part along the y direction. In the present example, four elements along 
the cantilever are found to be sufficient to describe heterogeneous effects, 
and the projection is made independently for each line across the beam. 
The figure El shows the resulting displacement field, which clearly exhibits 
an heterogeneous behavior across the cantilever. One should emphasize, 
that the free end of the cantilever experiences a phase change of 0.8 rad 
(Figll3p. This value is to be compared to the one arising only from a 12 nm 
height modification (Fig. 1140 . which is 0.18 rad (first term in Eq.()51|)) ; thus 
demonstrating the necessity of taking both height and slope changes into 
account. Moreover, the displacement field shows the cantilever part close to 
the anchoring rises up as the temperature increases, whereas the cantilever as 
a whole bends down. This is thought to be the signature of an under-etched 
cantilever, as the remaining silica pedestal dilates and pushes the surface 
up. This example typically illustrates the extensive amount of informations 
provided by full-field measurement compared to pointwise ones. 

5 Conclusions 

A Nomarski imaging interferometer is used to measure the differential topog- 
raphy of reflective objects. A phase modulation is introduced to measure the 
phase map arising from the object with a sensitivity independent on the ac- 
tual phase value. The phase measurement is shown to be shot-noise limited, 
and a measurement reproducibility of almost 10 pm is achieved. The use of 
four integrating buckets with a least-square algorithm is shown to improve 
the phase retrieval in case of low fringe constrast. Moreover, this allows to 
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reconsider the phase unwrapping problem, so that it is easy to deal with 
highly curved objects. 

The physical origin of the measured optical phase is calculated, and ex- 
hibits that the phase difference is due to both height and slope variations. 
A simple calibration procedure allows one to check for the phase sensitiv- 
ity to slope variations. These two effects are then decoupled expanding the 
displacement field onto a chosen functions basis. 

The out-of-plane displacement field of a MEMS cantilever subjected to 
bimaterial effect is then recovered thereby proving the ability of such a set-up 
to provide a reliable full-field kinematic measurement without surface mod- 
ification. This tool provides then a full displacement field using a common- 
path interferometer. Moreover, its stability, as well as its ability to operate 
through a wide range of media make this set-up a very powerful tool for 
the study of the mechanical behavior of micro-electro-mechanical systems, 
providing an extensive amount of reliable information. 
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A Appendix : Effect of the numerical aper- 
ture 



The section 22.41 is devoted to the calculation of the optical phase difference 
arising from the topography of the object under scope. This difference de- 
pends on the incidence of the ray on the surface. When using an imaging 
system, this difference depends therefore on the numerical aperture of the 
objective lens. The expression ()18|) was obtained under the assumption that 
the incidence angle is low enough. To assess this assumption's validity, let us 
consider the contribution of each incidence in the figure's plane to the total 
intensity in Eq0] 



dl(a) = ( dlo + gL4cos 



4nnAz 

— cos(a) + if) 

A 



da (60) 



where if) stands for the optical phase independent from the topography. As 
the fringe spacing depends on the angle a between the ray and the optic axis, 
the optical phase arising from the step in Fig. |5] is obtained by weighting 
and summing the contributions of each rays impinging on the sample at a 
point of the field of view 



I + A—^- / cos 

sin {a max ) J 



4nnAz 

cos(a) + ip 



P(a) 2 sm(a)da (61) 



A 

otmax is related to the numerical aperture of the objective lens NA 

NA = n sin(a max ) (62) 

and to the apodization function P(a). The choice of this function has been 
widely discussed |T[JJ 1201 EE] and let us the assume 

P(a) = (cos(a)) m (63) 

where m is a parameter used to describe the apodization effect. The equation 
(JU) turns into 

/ = J + A m F NA , m (Az, ^j) (64) 
for m = (Herschel condition), 

^ / „ ,s sin(kAz(l — cos(a max ))) ., . . . .-,,„■, 

F NAfi (Az^) ~ \ ^ V 7^ 008 ^1+008 0^ +^ 65 

kAz(l - cos(a maa; )) 
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where k = -y-. If m ^ 0, one can derive the Taylor expansion of ()64j) with 
respect to a max 

FNAm(Az,ijj) = — 27 7 / cos(2A;A2; cos(a) + if))P{a) 2m sin (a) da 

sm {a max ) J 

( 1 — 2m 1 
K cos(2^ + «) + (_ co S (2*A. + ,) + - S in(2*Az + ,)*A, 

f 1 + A(m 2 - m - k 2 Az 2 ) . , . 
+ f i — ^ cos(2A;Az + ^) 

kAzil - 4m) . . ,.\ 4 
+ ^2 sm ( 2A;A2 + ) «L + • • • 

Let us define the fringe spacing iinterf, as the value of kAz where the intensity 
(J64j) reaches its first local maximum. As the numerical aperture increases, 
iinterf 7^ 7T- Let us thus define the gap 

tinterf = 1 — ~ (67) 

Computing e inter f, when the numerical aperture ranges from to 0.6, and m 
ranges from to 2 shows that e inter f < so that the fringe spacing increases 
with the numerical aperture (see equation ([67)1 ). if the numerical aperture 
is less than 0.4, €i nte rf no longer depends significantly on m, so that the 
numerical aperture value is sufficient to retrieve Az from the intensity value. 
The above correction of the equation (JTBJ) is easy as long as FNA,m(Az,ip) is 
pseudo-periodic in a Az range sufficient to describe the topography. This is 
assessed by defining the ratio 

(68) 



where iw is the kAz value for which Fz\M, )m (Az, if)) reaches its tenth local 
maximum If the numerical aperture is less than 0.3, is 1, thereby proving 
that F N A,m(Az, if)) is pseudo-periodic on the defined range. The correction 
of Eq. (fTHj) then reads 

AimAz , 
0>< = — <«" 
where l depends on the numerical aperture and the apodization function. 
The figure shows the value of t as a function of both the numerical aperture 
and m. If the numerical aperture is less than 0.3, i no longer depends on m, 
and therefore allows for a direct retrieval of the phase (f>h- 
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Figure 2: Typical interference pattern obtained in water with two 70 x 20 x 
0.84/im 3 microcantilevers and a shear distance d ~ 50/im (NA=0.3). 
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Figure 3: Schematic view of a Wollaston prism. 
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Figure 4: Ray tracing for a plane object. 
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Figure 5: Ray tracing in the case of a tilted and stepped sample (height 
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Figure 6: Ray tracing for a sample subjected to slope variations. 
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Figure 7: Calibration measurement of the mean phase induced by a sample 
function of its tilt. 
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Figure 8: Typical wrapped phase map obtained in water with two 70 x 20 x 
0.84/xm 3 microcantilevers and a shear distance d = 53.4/im (NA=0.3). 
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Figure 10: Typical unwrapped phase map obtained in water with two 70 x 
20 x 0.84/im 3 microcantilevers and a shear distance d = 53.4/im (NA=0.3). 
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Figure 11: Typical experimental phase noise probability density (dots) and 
its fit by a Gaussian distribution (solid line). 
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Figure 12: Estimation of the reproducibility on the measurement of a differ- 
ential topography as a function of the exposure time. 
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Figure 13: Measured phase map change when the cantilever is subjected to 
bimaterial effect. 
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Figure 14: Displacement field calculated from the measured phase map 
change shown in FigfEfl 
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Figure 15: Evolution of the correction factor l when the numerical aperture 
ranges from to 0.4, and the exponent m ranges from to 2. 
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